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ABSTRACT 

This paper presents a detailed theoretical study of the two-point correlation 
function ^ for both dark matter halos and the matter density field in hve cosmological 
models with varying matter density Dm and neutrino fraction The objectives 
of this systematic study are to evaluate the nonlinear gravitational effects on C to 
contrast the behavior of ^ for halos vs. matter, and to quantify the redshift evolution 
of ^ and its dependence on cosmological parameters. Overall, ^ for halos exhibits 
markedly slower evolution than ^ for matter, and its redshift dependence is much 
more intricate than the single power-law parameterization used in the literature. Of 
particular interest is that the redshift evolution of the halo-halo correlation length 
ro depends strongly on Dm and D^/, being slower in models with lower Dm or higher 
Djy. Measurements of ^ to higher redshifts can therefore be a potential discriminator 
of cosmological parameters. The evolution rate of ro for halos within a given model 
increases with time, passing the phase of hxed comoving clustering at z ~ 1 to 3 
toward the regime of stable clustering at z ~ 0. The shape of the halo-halo on the 
other hand, is well approximated by a power law with slope —1.8 in all models and is 
not a sensitive model discriminator. 
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1. Introduction 

The two-point correlation function ^(r) is a fundamental quantity in cosmology. For the 
density field of matter, it provides the most basic statistical measure of gravitational clustering, 
and its Fourier pair is simply the power spectrum P{k). For galaxies, ^ offers a convenient measure 
of their clustering strength, and a determination of ^ is among the most important goals of every 
major galaxy survey. Although ^ is not independent from P, the estimators for the two quantities 
have different error properties and therefore both carry useful information. 

The two-point correlation function of nearby galaxies in physical space has been determined in 
several galaxy redshift surveys and is found to be well approximated by a power law, ^ = (r/rg 
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over the distance range 1 < r < 20/i“^ Mpc. Analysis of the early CfA sample of 2400 galaxies 
(complete to ms = 14.5) gives a correlation length of ro = 5.4 ± 0.3Mpc and a slope of 
7 = 1.77 for 1 ^ r ^ 10Mpc (Davis & Peebles 1983). Subsequent optical surveys generally 
have yielded consistent results. The combined CfA2 and Southern Sky Redshift Surveys of 12812 
galaxies, for example, give ro = 4.95 ± 0.13Mpc and 7 = 1.73 ± 0.03 when cluster galaxies 
are removed (Marzke et al. 1995). The Las Campanas Redshift Survey (Shectman et al. 1996) of 
19558 galaxies is found to have ro = 5.06 ± 0.12Mpc and 7 = 1.862 ± 0.034 (Jing, Mo, and 
Borner 1998). Surveys at other wavelengths have shown larger variations. The 1.2 Jy IRAS survey 
of 10000 galaxies yields ro = 3.76Mpc and 7 = 1.66 for r < 20Mpc (Fisher et al. 1994). 
A small sample of 300 1.4 GHz radio sources is found to have ro = 11 h~^ Mpc and 7 = 1.8, where 
the enhanced correlations are attributed to the high optical luminosity and rich environments 
inhabited by radio galaxies (Peacock & Nicholson 1991). 

Although measurements of ^(r) for nearby galaxies are abundant, the two-point correlation 
function at higher redshifts is an important quantity that has only begun to be determined. Deep 
redshift surveys are clearly required for this purpose, and initial results have been somewhat 
conflicting. The Autofib survey of 1100 galaxies from 33 pencil beams, for example, shows 
no evidence of evolution, giving ro = 6.2 ± 0.4Mpc and rg = 6.8 ± 0.4Mpc for the 
samples below and above the survey median redshift {z) = 0.16, respectively (Cole et al. 1994). 
The Canada-France Redshift Survey, on the other hand, finds a decreasing correlation function 
with increasing redshift, with ro = 1.33 ± 0.09Mpc at z = 0.53 (Le Fevre et al. 1996). A 
sample of 183 field galaxies covering 216 square arcmin from the CNOC (Canadian Network for 
Observational Cosmology) cluster survey yields ro = 2.2 ± 0.5 h~^ Mpc for (z) = 0.37 (Shepherd 
et al. 1997). A sample of 248 K-band selected galaxies from two fields of combined area of 27 
square arcmin gives ro = 2.9 ± 0.34 h~^ Mpc at z = 0.34, ro = 2.0 ± 0.3 h~^ Mpc at z = 0.62, 
and ro = 1.4 ± 0.21 h~^ Mpc at z = 0.97 (Carlberg et al. 1997). The shape of ^ shows little 
variation with redshift and is well fit by a power law with a slope 7 = 1.7 to 1.8 in all surveys. 
More measurements from deep surveys covering larger areas of the sky will provide more robust 
determination of ^(r) beyond z ^ 0. 

On the theoretical side, models for the redshift evolution of ^(r) have been simplistic for the 
most part, and the interpretation of observational results is often cast either in the regime of 
stable clustering in physical space (^ oc a^, where a = (l-|-z)“^), or fixed clustering in comoving 
space. Uncertainties and statistical fluctuations in older data sets due to small galaxy numbers 
and limited survey fields are often too large to allow a finer distinction beyond these two regimes, 
but data from the ongoing and planned future deep redshift surveys certainly have the potential 
for discriminating among cosmological parameters, if the evolution of ^(r) follows distinct patterns 
in different models. This provides the motivation for the extensive investigation of ^ presented 
here for a wide range of redshift and cosmological models. 

Theoretical predictions of ^ fall in two categories; the density field and dark matter halos. 
Much work has been done for the former (Jain 1997; Sheth and Jain 1997; Matarrese et al. 1997; 
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Jenkins et al. 1998; and references therein), particularly in the context of scale-free models with 
power-law P{k) and the cold dark matter (CDM) models. No systematic study, however, has 
been carried out to investigate a spectrum of models with different matter densities as well as 
neutrino masses over a large redshift range. Theoretical calculations of ^ for the density field can 
enhance our general understanding of gravitational clustering, but the results cannot be compared 
directly with the observed ^ without the knowledge of the relationship between the density field 
and galaxies. 

The two-point correlation function of dark matter halos, on the other hand, bears closer 
resemblance to the observed Analytical bias models have been constructed to relate ^ for 
the density field to halos (Moscardini et al. 1998), but a more direct, less model-dependent 
approach is to compute the halo-halo ^ from numerical simulations of structure formation. Since 
this requires more extensive effort, only a few individual models have been studied. Gelb & 
Bertschinger (1994b), Brainerd & Villumsen (1994), Bagla (1997), and Jing (1998), for example, 
examined ^ of both matter and halos in the standard CDM model (although with different halo 
identification schemes. See § 3.2). The correlation function in the = 0.3 coldj-hot dark matter 
(C-I-HDM) model has been studied in Klypin et al. (1993) and Jing et al. (1994), but this model 
has been shown to possess too little small-scale power at high redshifts to account for the amount 
of observed collapsed objects. Colin et al. (1997) is the only work thus far that has specifically 
quantified the redshift evolution of halo-halo ^ in different models. Both Dm = 1 and Dm = 0.2 
CDM models were examined, but their choice of high Hubble constant, h = 1, exceeds all current 
observational values and also gives the Dm = 1 CDM power spectrum a shape parameter of T = 1, 
which is far above T ~ 0.25 indicated by galaxy clustering and the value of T = 0.5 for the 
standard CDM model. 

In this paper I present results from an extensive theoretical study of the properties and 
evolution of ^{r) in three classes of structure formation models that are of current cosmological 
interest: CDM, C-I-HDM, and fiat low-density CDM models with a cosmological constant (LCDM). 
The objectives are to (1) evaluate the nonlinear gravitational effects on (2) compare the behavior 
of ^ for matter vs. halos, which reffects the bias between matter and halos and is sometimes a 
confused issue in the literature; and (3) quantify the redshift evolution of ^ and its dependence on 
cosmological parameters. No such systematic study over this wide range of cosmological models 
and redshift has been carried out thus far. In order to determine the nonlinear evolution of 
numerical simulations of structure formation in five flat models with matter density Dm = 1, 0.5, 
0.3, and neutrino fraction D^, = O.I and 0.2 are performed (see § 2). Results of ^ for both the 
density field and halos from linear as well as nonlinear calculations are discussed in § 3. Particular 
emphasis is placed on quantifying and contrasting the redshift dependence of the halo-halo ^ in 
different cosmogonies, which can serve as a potential discriminator of cosmological parameters 
when more high-redshift data are collected. 
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2. Simulations of Nonlinear Clustering 

Numerical simulations are performed to determine the nonlinear two-point correlation 
functions for both the density field and dark matter halos in five cosmological models that include 
two C-I-HDM, two LCDM, and the standard CDM model. For this study to form a controlled 
numerical experiment, the simulations for the C-I-HDM and LCDM models reported here are 
all performed in a (100 Mpc)^ comoving box with a Plummer force softening length of 50 kpc 
comoving (for the given value of the Hubble constant in each model), and identical phases are used 
in the initial conditions for the runs. Two additional large simulations in a (640 Mpc)^ box with a 
force softening length of 160 kpc and a different set of initial conditions are used to test the effects 
of numerical resolution and realization on our results (see § 3). The primordial power spectra for 
all models have a spectral index of n = 1, and the density fluctuations are drawn from a random 
Gaussian distribution. The gravitational forces are computed with a particle-particle particle-mesh 
(P^M) code (Bertschinger & Gelb 1991; Gelb & Bertschinger 1994ab; Ma & Bertschinger 1994b). 
The number of simulation particles used to represent the cold dark matter is 256^ for the 640 Mpc 
CDM run, and 128^ for the rest. For the hot dark matter component in C-I-HDM models, 128^ 
and 10 X 128^ particles are used in the = 0.1 and 0.2 models, respectively. Although the latter 
is needed for a fine sampling of the neutrino velocity phase space (Ma &: Bertschinger 1994a), the 
two-point correlation function - the subject of study here - mostly probes the spatial distribution 
of matter. Tests with a reduced number of hot particles have been performed in the Q.i, = 0.2 
model to verify that 128^ is sufficient for an accurate determination of 

For the C-I-HDM models, only = 0.1 and = 0.2 are studied here because structure 
forms too late in models with Hi, > 0.2 to account for high-redshift galactic structures (Ma et ah 
1997 and references therein). A baryon fraction of Hb = 0.05 and a Hubble constant of /i = 0.5 
are assumed, and the cold dark matter fraction is Hcdm = 0.85 and 0.75 for the two models, 
respectively. The rms linear mass fluctuation in spheres of radius 8 h~^ Mpc is a% = 0.9 and 0.82 
for Hi, = 0.1 and 0.2, when the 4-year COBE result (Bennett et al. 1996; Gorski et al. 1996) is 
used as the normalization. 

The two LCDM models chosen for the simulations have (Hm,HA,h) = (0.3,0.7,0.75) and 
(0.5,0.5,0.7). The parameter F = Hm h characterizing the shape of the power spectrum is 0.23 and 
0.35 for the two models. Both models are COBE-normalized, with erg = 1-28 and 1.53, respectively. 
The COBE-normalized standard CDM model with Hm = 1 and h = 0.5 (no baryons) is also 
included in this study as a point of standard reference despite of its excess power on small scales. 
The corresponding fTg = 1.4 is a factor of ~ 2 higher than that required by cluster abundance, but 
the results quoted here can be easily translated to a lower normalization by rescaling the redshift. 
The results for the standard CDM model are computed from a P^M simulation performed by Gelb 
& Bertschinger (1994ab), which employed a numerical algorithm and numerical parameters similar 
to the LCDM and C-I-HDM runs used here; 66 kpc for the Plummer force softening distance, 100 
Mpc box, and 144^ particles. The baryon fraction is set to zero for both LCDM and the standard 
CDM models so that the effects of changing Hm can be clearly studied. 
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3. Two-Point Correlation Function 

The two-point correlation functions ^ for the mass density field and dark matter halos have 
different shapes and follow distinct evolutionary histories. They are discussed separately in § 3.1 
and § 3.2. Results for the evolution of the correlation length are discussed in § 3.3. All figures for 
^ are presented in physical coordinates r, but its behavior in comoving space is also discussed. 


3.1. ^ of the Density Field 

The two-point correlation function of the density field measures the lowest-order gravitational 
clustering of matter in a given cosmological model. It is the Fourier transform of the power 
spectrum P{k) for the density fluctuation 5p/p, and can be written as 

«r)= (1) 

JO K kr 

In the linear regime, ^ is easily calculated from the linear P{k), which can be obtained from either 
direct integration of the coupled Boltzmann and Einstein equations, or analytical approximations 
for various cosmological models that have been published in the literature (e.g., Bardeen et al. 
1986; Efstathiou et al. 1992; Ma 1996). In the nonlinear regime, ^ can be calculated either from 
the particle positions in numerical simulations directly, or from integration of the few published 
fitting functions for the nonlinear P{k). Eor the nonlinear P{k), Hamilton et al. (1991) examined 
scale-invariant models with power-law spectra, Jain et al. (1995) and Peacock & Dodds (1996) 
studied the flat and open CDM models, while Ma (1998) provided a general approximation valid 
for CDM, LCDM, as well as C-I-HDM models. 

Figures 1 to 3 compare ^(r) for the linear and nonlinear density fields at various redshifts for 
the five models chosen for this study - the standard CDM (Fig. 1), the Dm = 0.5 and 0.3 LCDM 
models (Fig. 2), and the = 0.1 and 0.2 C-I-HDM models (Fig. 3). The dashed curves show the 
linear ^ computed from equation (|^). For the CDM and LCDM models, the input linear P{k) is 
the Bardeen et al. fitting function (1986) with the shape parameter F = VL^h (no baryons are 
assumed in these models). For the C-I-HDM models, the input is the fitting function of Ma (1996) 
for the density-averaged linear P{k), where P{k) = -|- (1 — D^,) measures the total 

gravitational fluctuations contributed by the separate hot and cold components. 

The circles in Figures 1 to 3 show the nonlinear ^ calculated directly from the particles in 
the (100 Mpc)^ simulations described in § 2. As a test of the robustness of our numerical results, 
the solid squares in the bottom panels of Figures 1 and 2 compare ^ computed from different 
simulations of the CDM and the Dm = 0.3 LCDM models. Both simulations are performed in a 
much larger box - (640 Mpc)^ in comoving volume - and a force softening distance of 160 kpc 
(i.e., 80/i“^ kpc for CDM and 120kpc for LCDM). The CDM simulation uses 16.8 million 
particles while the LCDM has 2.1 million. As expected, the squares fall slightly below the circles 
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at small separation due to the larger force softening in the big-box simulations, whereas the circles 
fall below the squares at comoving separation x 10 h~^ Mpc due to their smaller box size. In the 
intermediate regime, however, the results agree well between the two runs that have very different 
particle numbers, box sizes, and force resolution. This demonstrates that our determination of ^ 
is robust within the resolution range of the simulations. 

As discussed earlier, the nonlinear ^ can also be computed from equation (|^ if the 
corresponding nonlinear power spectrum is known. The solid curves in Figures 1 to 3 show the 
results for all five models when the analytical approximation of Ma (1998) for the nonlinear P{k) 
is used as input in equation (||). The results agree well with the nonlinear ^ computed directly 
from the simulation particles, providing a consistency check. 

As Figures 1-3 illustrate, the effects of nonlinear gravitational clustering drastically change 
the shape of the linear ^ below about comoving Mpc. At separations of ~ 0.1comoving 

Mpc, the nonlinear contribution boosts the linear correlation amplitude by more than one order 
of magnitude. As a result, the nonlinear effect on the shape of ^ is to straighten up the linear 
curves partially, but as the solid curves show, the nonlinear ^(r) is far from being a simple power 
law. In particular, ^ acquires a steeper slope in the intermediate scale around 1 h~^ comoving 
Mpc in all five models. Such a change of slope is consistent with similar behavior in the nonlinear 
power spectrum in Hamilton et ah (1991), Jain et al. (1995), and Ma (1998). It is also reported in 
scale-free and CDM models by Padmanabhan (1996), Jain (1997), and Jenkins et al. (1998). This 
non-power-law behavior of ^ of the density field, however, should not be viewed as an inconsistency 
with observations. The ^ computed from halos discussed in the next section bears closer reality 
to the observed ^ for galaxies, and as will be shown, the shape of halo-halo ^ is indeed well 
approximated by a power law. 

In the intermediate range between 1 and 10 comoving Mpc, Figures 1 to 3 show the 
interesting feature that the nonlinear ^ dips below the linear ^ before rising rapidly above it on 
smaller scales. The decrement is present only at z ^ 1, but it becomes more prominent and 
extends over a larger range of r as time goes on. This feature is not due to numerical artifacts 
in the simulations (such as small box size or statistical fluke) because it is seen in all models, in 
both 100 Mpc and 640 Mpc simulations with different resolution and realization, and in both ^ 
computed from the particles (i.e. circles and squares) and that computed from integration of the 
analytical nonlinear P{k) (i.e. solid curves). Mathematically, despite the fact that the nonlinear 
P{k) is always larger than or equal to the linear P{k) (within the < 10% fitting error of Ma 
(1998)), the nonlinear ^ can fall below the linear ^ for a limited range of r because the factor 
sm{kr)/kr in equation (|T]) is oscillatory. Physically, the decrement may reflect the depletion of 
matter on quasi-linear scales that is fueling the process of gravitational infalls in more clustered 
regions on smaller scales. The decrement also implies that the nonlinear correlation lengths evolve 
more slowly than the linear predictions, which contributes to the slower evolution of the halo-halo 
^ to be discussed in the next section. This interesting feature may warrant further analytical 
study. 
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The dependence of ^{r) on the expansion factor a is related to the growth rate of Pik), a well 
understood phenomenon in the linear regime. For the standard CDM model, the linear P and ^ 
both grow as a^. For C+HDM models, the linear P and ^ also increase as above the neutrino 
free-streaming scale, but the growth is retarded below this scale because the phase-mixing of 
neutrinos reduces the density perturbations. An accurate analytic approximation (with < 10% 
error) for this scale-dependent effect is given by equation (11) of Ma (1996). For critically-flat 
LCDM models, the linear P and ^ again grow as a? until the Universe becomes A-dominated 
at 1 + z K, Um ' . After this epoch, the growth is slowed down (Heath 1977). However, unlike 
the C-I-HDM models in which the neutrino free-streaming imprints a characteristic scale, the 
retardation of the linear growth rate in the LCDM model is independent of length scales. In the 
nonlinear regime, the evolution of ^ for the density field in the CDM model has been discussed in 
Jain (1997) and Matarrese et al. (1997). 


3.2. Halo-Halo ^{r) 

In order to determine the correlations of the simulated halos, the Denmax algorithm (Gelb 
& Bertschinger 1994ab; Ma & Bertschinger 1994b; Frederic 1995ab) is first used to identify dark 
matter halos in the output of the simulations. The Denmax algorithm is designed to associate 
particles with density peaks by moving the particles along the gradient of the mass density field 
until they settle into the peaks. The density field is first computed on a very fine grid from the 
particle positions using the triangular-shaped cloud interpolation scheme (Hockney & Eastwood 
1988), and then convolved with a Gaussian filter to remove grid effects. After the particles are 
associated with density peaks by Denmax, the energy of the particles in each halo is computed 
and the unbound ones are removed. An overdensity of 200 is also imposed to select virialized 
halos that are potential galaxy hosts. Due to the small smoothing length required for identifying 
galactic-size halos (typically 1/1000th the box length), this method is more CPU-intensive than 
the friends-of-friends (FOF) algorithm used in previous numerical studies of halo-halo correlation 
functions (e.g., Brainerd & Villumsen 1994; Colin et al. 1997). The Denmax algorithm is chosen 
here, however, because extensive tests have demonstrated that it generally identifies distinct 
density concentrations more reliably than the FOF algorithm, for which an arbitrary linking 
length must be assigned (Gelb & Bertschinger 1994ab; Frederic 1995ab). 

All dissipationless simulations suffer to some degree the overmerging problem in regions of 
high overdensities. When individual dark matter halos merge in these simulations, they form a 
larger, smooth halo with little substructure due to the lack of dissipation. The extent of this 
problem can be quantified by examining output of hydrodynamic simulations, which are designed 
to follow both the dark matter and a gas component that is allowed to cool and condense to high 
densities. Massive dark matter halos of the size of galaxy clusters in these simulations have been 
seen to host multiple subclumps of condensed gas - presumably the sites of galaxies - that had 
survived the merger due to their high densities and small cross sections (e.g., Katz et al. 1992; 



Evrard et al. 1994). It is not easy to identify galactic halos in such highly clustered regions from 
the dark matter distribution alone, and various methods have been proposed to locate galactic 
sites in clusters in iV-body simulations. However, tests against hydrodynamical simulations have 
so far indicated that no method proves satisfactory on cluster scales (Summers et al. 1995). 

Outside of the dense cluster environment, on the other hand, dark matter halos identified 
by reliable algorithms in iV-body simulations have been found to give robust predictions for 
galaxy clustering statistics. For the two-point correlation function, for example, Summers et al. 
(1995) found good agreement between ^ determined from dark matter halos and from the gas in 
hydrodynamical simulations on scales above 1 Mpc, and most halo identification methods gave 
consistent results (see their Fig. 12). Gelb & Bertschinger (1994b) studied two break-up methods 
that were designed to reconstruct the galactic substructure in clusters by splitting the massive 
halos in A^-body simulations into multiple sub-halos. The amplitude of the two-point correlation 
was indeed boosted significantly on scales below 1 Mpc due to the break-up procedure. On scales 
above 1 Mpc, however, both the shape and the correlation length were hardly changed whether 
the massive halos were split up or not (see their Figs. 9 and 10). These studies all indicate that 
the overmerging problem in A^-body simulations is confined to scales of ~ 1 comoving Mpc and 
below; it therefore does not affect the determination of the correlation length (discussed in the 
next section). Moreover, most measurements of ^ are based on surveys of galaxies in blank fields 
devoid of massive clusters (e.g. Shepherd et al. 1997; Le Fevre et al. 1996; Carlberg et al. 1997). 
The problem with galaxy identification in the cluster environment therefore should not be of 
serious concern to this study. 

Figures 4-6 show the resulting two-point correlation functions of the resolved halos determined 
from our simulations. Five models - the standard CDM (Fig. 4), the = 0.3 and 0.5 LCDM 
models (Fig. 5), and the = 0.1 and 0.2 C-I-HDM (Fig. 6) - are shown for comparison. For 
each model, the figures illustrate the redshift evolution by plotting ^ from z = 0 up to 4. The 
left and right panels are computed from virialized halos in two mass ranges, M > 10^^'^ Mq 
and M > 10^^'^ Mq, respectively. These two different mass limits are chosen to illustrate the 
dependence of ^ on halo populations. The figures show that the more massive halos have a higher 
correlation amplitude and length. Since the observed ^ is determined from field galaxies not 
associated with massive clusters (as discussed above), I have tested the effects of massive halos 
on the determination of ^ by imposing an upper mass cut of M < 10^^ Mq. Due to the rarity of 
these high-mass objects, ^ is found to depend little on whether these objects are included in the 
sample or not. For simplicity, such an upper mass cut will therefore not be imposed in the results 
reported here. 

As Figures 4-6 show, the shapes of the halo-halo ^ do not depend noticeably on the models 
or the halo mass cuts. In contrary to ^ for the density field (Figs 1 to 3), the halo-halo ^ is 
well approximated by a power law with the observed slope of ~ —1.8 between 1 and 10 
comoving Mpc. On scales of below ~ 1 h~^ comoving Mpc, starts to flatten out as z approaches 
0 in the CDM and Dm = 0.5 LCDM models that have relatively large small-scale power. This 



- 9 - 


behavior largely reflects the over-merging problem in dissipationless simulations discussed above. 
The flattening of ^ disappears when either halo break-up procedures are applied (e.g., Gelb & 
Bertschinger 1994b) or gas particles are used as galaxy tracers. Since the shape of ^ carries little 
information about the underlying cosmogony, no attempts will be made here to correct for the 
flattening of ^ on small scales. Instead, the focus of this study is on the more interesting quantity, 
the correlation length and its evolution. They are determined by the reliable part of ^ at r > 1 h~^ 
Mpc and are largely unaffected by the overmerging problem. 


3.3. Redshift Evolution of ^ and vq 

Comparison of the two-point correlation function for halos in Figures 4-6 to that for the 
density held in Figures 1-3 shows that the evolution of the halo-halo ^ is markedly slower than the 
matter density held in all models. Relative to the density held, the clustering amplitude of halos 
is higher at z ~ 4, but it invariably drops below that of the density field by z = 0 due to the slow 
evolution. Such “anti-biasing” effect has been noted in previous studies of the CDM model (Gelb 
& Bertschinger 1994b; Colin et al. 1997); the present study shows that this effect occurs in LCDM 
and C-I-HDM models as well. It is attributed to the continuous disruption and mergers of halos at 
late times that reduce the clustering of halos to below that of the mass. The relationship between 
the density field and the halos shows complicated dependence on both length and time scales and 
on cosmological parameters. Care must therefore be taken not to compare ^ for the density field 
to the observed ^ since the former does not reflect directly the observed clustering properties of 
galaxies. 

The evolution of ^ for halos depends strongly on the cosmological parameters. Figure 7 
shows the redshift dependence of the physical correlation length rQ{z) for the simulated halos for 
two mass ranges in all five models. The fastest growth occurs in the pure CDM model, in which 
the physical correlation length ro increases by a factor of 8 from z = 4 to 0. The correlation 
length changes more mildly in the LCDM and C-I-HDM models. For instance, ro increases by 
a factor of 3 from z = 3 to 0 in the = 0.2 model and a factor of 5.4 from 2 : = 4 to 0 in 
the Dm = 0.3 LCDM model. When expressed in comoving coordinates, the halo-halo ^ in fact 
decreases with time at higher redshifts and is seen to rise only at a later time (see also Brainerd 
and Villumsen 1994; Bagla 1997). However, one should be careful not to interpret this behavior 
as actual decrease in the physical clustering amplitude at early times: as Figures 4 to 7 clearly 
illustrate, the physical ro and ^(r) always increase monotonically with time in all models. Instead, 
the decrease merely reflects the fact that the halo clustering rate is slower than the state of fixed 
clustering in comoving coordinates early on. It eventually surpasses fixed comoving clustering 
when the comoving correlation length begins to increase (see Fig. 8 below). 

The top panels of Figure 7 illustrate that ro increases more slowly in models with smaller 
Dm. This feature is consistent with the slower growth of structure in low-density models. The 
lower panels compare different C-I-HDM models and indicate that larger D^^ models show slower 
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evolution. This can be explained by the slower growth due to the large neutrino free-streaming 
effects in high-fl^y models. The left and right panels of Figure 7 illustrate that vq also depends 
on the minimal halo mass selected for the sample, and is larger for the more massive halos (right 
panels) due to their high clustering amplitudes. 

To quantify the evolution of the two-point correlation function with redshift, I find it most 
convenient to parameterize the physical and comoving correlation lengths directly as 


roiz) = (l + z)-'^(^)ro(0), 

xo{z) = (l + z)-'’(^)+ixo(0). (2) 

The parameter rj itself increases with time because as Figure 7 illustrates, the slope of ro(z) 
generally steepens as z decreases, and the rate of change depends on and The form of 
equation (| 2 |) is chosen to improve the power-law parameterization commonly used in the literature, 
where both the shape and redshift dependence of ^ are written as power laws with exponents 7 
and e: 


C{r,z) 

C(r,z) 



(l + z)-^-^Cir, 0 ), 


( 3 ) 


and r is the pair separation in physical coordinates. It follows from equation (^) that the physical 
and comoving correlation length evolve as 

ro{z) = ( 1 z)^^ro(O), 

xo(z) = {l + z) xo(0). (4) 

As Figure 7 illustrates, however, the power-law redshift dependence in equation (Q) does not hold 
over the redshift range 0 to 4 for any model. Moreover, since the parameter e in equation (^) 
was introduced to characterize the redshift evolution of ^ and not rg, the evolution of tq depends 
inconveniently on both e and 7 . The simple form of equation (||) will therefore be used here. To 
translate r] into the parameter e used in the literature, use 


e = j rj — 3 . 


( 5 ) 


Some useful theoretical predictions for the evolution of ^ are: (1) The linear theory predicts 
^ oc (1 -|- z)~'^ for the standard CDM model, corresponding to e = 7 — 1, or e ~ 0.8 and rj ~ 2.1 
for 7 w 1 . 8 ; ( 2 ) clustering fixed in comoving coordinates corresponds to xq{z) = xo(0), implying 
r/ = 1, or e w —1-2; (3) stable clustering in physical space corresponds to ^ oc (1 -|- z)~^, and 
therefore e = 0, or r/ ~ 1.7. 

Figure 8 shows r]{z) for the simulated halos in our study. The rise of r]{z) with time indicates 
that the halo correlation length rg grows faster than a simple power law with a constant e assumed 
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in the literature (as seen in Figure 7). The parameter rj{z) also depends on both cosmological 
parameters and The label on the right along the vertical axis marks the corresponding 
e when 7 = 1.8 is assumed for the slope of When r/ crosses unity in a given model, it 
indicates that gravitational clustering has become faster than no evolution in comoving space (i.e., 
xo{z) = xo(0)). This epoch corresponds to the turn-around redshift at which the initially falling 
^ first starts to increase in comoving space. At 2; ~ 0, the CDM, flm = 0.5 LCDM, and fli, = 0.1 
C-I-HDM models all approach the stable clustering regime of t/ ~ 1.7, whereas the = 0.3 
LCDM model exhibits less rapid evolution of t/ ~ 1.3. Stable clustering has also been seen for 
the density field on small scales in the pure CDM model and models with scale-free spectra (Jain 
1997). Our study here shows that halos also approach stable clustering in high-Dm models. 

The only detailed numerical study of the evolution of halo-halo ^ thus far is Cohn et ah 
(1997), in which the Dm = 1 and Dm = 0.2 CDM models are examined. But as discussed in §1, 
the Dm = 1 CDM results reported here cannot be compared directly to theirs because of their 
choice of high Hubble constant h = 1 and large value of the shape parameter T = 1. Since the 
standard CDM model with T = 0.5 used in this paper already suffers from the well-known excess 
small-scale power, a model with T = 1 would produce even more power and worsen the problem. 
This may be the main cause of the steep slope 7 ~ 2.2 to 2.4 for ^ in Table 2 of Colin et ah, 
which is significantly larger than the value of 7 ~ 1.8 found here and in other works. One should 
also note that because ro{z) depends on both 7 and e (see eq. @), their high values of 7 lead to 
the large positive e in their Figure 9 (where e ranges from ~0atz = 5 tol. 5 at 2 : = l). When 
expressed in terms of the parameter rj introduced earlier, their large values of 7 and e do cancel 
out to some extent, giving rj ^ 1.3 at z = 5 and 1.87 at 2 = 0, but it is still significantly larger 
than our e in Figure 8 . 


4. Conclusions and Discussion 

Motivated by the extensive observational effort toward determining the evolution of the 
two-point correlation function ^(r), I have presented a detailed theoretical study of the the redshift 
dependence of ^ in five cosmological models with varying mass density Dm and neutrino fraction 
Dj,. The evolution is found to be much more intricate than the simple power-law parameterization 
of e commonly assumed in the literature, and a redshift-dependent variable rj is introduced to 
quantify the evolution of the correlation length: ro{z) = (1 -|- ?'o(0). The value of rj{z) 

depends strongly on Dm and D,y, and is smaller in models with lower Dm or higher Djy (Fig. 8), 
indicating a slower growth in tq. Within a given model, rj increases with time, passing the phase 
of fixed comoving clustering (77 = 1) at z ~ 1 to 3 toward stable clustering (rj ^ 1.7) at z ^ 0. 

Although the two-point correlation function ^ provides a fundamental statistical measure of 
gravitational clustering, care must be taken in specifying the underlying quantity that is being 
measured and in including the possible nonlinear effects. This study has demonstrated that ^ for 
the matter density field and dark matter halos have different shapes and evolve with different 
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rates. While the local slope of ^ for the nonlinear density field changes with length scales (Figs. 1 
to 3), the shape of the halo-halo ^ is well approximated by a power law with the observed slope 
of ~ —1.8 independent of cosmological parameters (Figs. 4 to 6). It is therefore the redshift 
evolution of the amplitude of ^ that will reveal the most about cosmogonies. Since halos generally 
correspond to galaxy sites (except in dense cluster environment), ^ for halos in simulations is a 
reasonable quantity to compare with observations. The correlation function for the density field 
provides important insight for theoretical studies of gravitational clustering, but it bears little 
direct relevance to the observed clustering of galaxies. 

On the observational side, ongoing deep redshift surveys of galaxies that cover larger areas of 
the sky should soon yield more robust determination of ^ beyond z ~ 0. Photometric redshifts 
determined from multi-color surveys and the two-point angular correlation function also offer 
complimentary measurements of ^ (Connolly, Szalay, & Brunner 1998; Postman et al. 1998; and 
references therein). We are currently analyzing galaxy clustering properties in two redshift surveys 
conducted with the Norris Spectrograph on the Palomar 200-inch telescope (Small et al. 1998). 
The combined surveys have more than 1000 galaxies in 0.2 < 2 : < 0.5. The implications of the 
theoretical results presented here will be discussed in Small et al. (1998). 

Although the evolution of ^ has already been counted as yet another strike against llm = 1 
models (Carlberg et al. 1997), it must be borne in mind that the nature of high-redshift galaxies 
has only begun to be unravelled, and that little is known about their dark matter halos. Due to 
the large uncertainties in the current observational results and theoretical understanding in galaxy 
evolution, this paper has chosen to investigate the evolution of ^ using simulated dark matter 
halos within well defined mass ranges. It is nonetheless conceivable that measurements of ^ at 
higher redshifts are probing different galaxy populations that are not fair counterparts of the local 
galaxies. Based on future observations, one should be able to refine the theoretical predictions 
made in this paper by incorporating different dark matter halo populations at higher redshifts. 


The author thanks Jim Frederic and Todd Small for valuable discussions. Supercomputing 
time was provided by the National Scalable Cluster Project at the University of Pennsylvania and 
the National Center for Supercomputing Applications. 
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Fig. 1.— Evolution of the two-point correlation function ^(r) (in physical coordinates) for the 
matter density field in the standard CDM model with = 1 and h = 0.5. The top panel shows 
results at as = 0.3, 0.5, and 1.0 (from bottom up), and the bottom panel shows us = 0.7 separately 
for clarity. (The COBE normalization corresponds to ag = 1.4 at z = 0, while ag ~ 0.7 is required 
to match the observed abundance of z ~ 0 galaxy clusters.) The dashed curves show linear-theory 
predictions. The symbols represent the nonlinear ^ computed from particles in the simulations - 
the circles are from a simulation in a (100 Mpc)^ box, and the solid squares for the fTg = 0.7 output 
is from a (640 Mpc)^ simulation. The solid curves, which agree well with the symbols, represent 
the nonlinear ^ from integration of the analytic nonlinear P{k) in Ma (1998). The dotted straight 
line indicates a slope of —1.8 in this and the following figures. 
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Fig. 2.— Same as Fig. 1 but for two LCDM models; flm = 0.5, Qa = 0.5, and h = 0.7 (top); and 
11m = 0.3, = 0.7, and h = 0.75 (bottom). For each model, four redshifts are shown: z = 4, 

3, 1.5, and 0 (from bottom up). The solid squares compare ^ computed from a large (640 Mpc)^ 
simulation of the flm = 0.3 model. 
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Fig. 4.— Evolution of the two-point correlation function ^(r) for the simulated dark matter halos 
in the standard CDM model with fljn = 1 and h = 0.5. Five epochs are shown; erg = 0.3 
(dotted), 0.5 (dot-short-dashed), 0.7 (short-dashed), 1.0 (long-dashed), and 1.4 (solid). (The CODE 
normalization corresponds to fig = 1.4 at z = 0.) The left and right panels show ^ computed from 
virialized halos in the mass range M > 10^^'^ Mq and M > 10^^'^ Mq, respectively. The dotted 
straight line indicates the observed slope of —1.8. 
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Fig. 5.— Same as Fig. 4 but for the two LCDM models: = 0.5, Qa = 0.5, and h = 0.7 (top); 

and rijn = 0.3, Ha = 0.7, and h = 0.75 (bottom). Five epochs are shown; z = 4 (dotted), 3 
(dot-short-dashed), 1.5 (short-dashed), 0.5 (long-dashed), and 0 (solid). 
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Fig. 6.— Same as Fig. 4 but for the two C+HDM models ~ = 0.2 (top) and Oi, = 0.1 (bottom). 

For Qiy = 0.2, the redshifts shown are; z = 2.6 (dotted), 1.5 (dot-short-dashed), 1.0 (short-dashed), 
0.26 (long-dashed), and 0 (solid). For the Qi, = 0.1 model, the redshifts are: z = 3 (dotted), 2 
(dot-short-dashed), 1.5 (short-dashed), 0.43 (long-dashed), and 0 (solid). (The dotted curves are 
omitted from the right panels since estimates of ^ are very noisy due to the small number of massive 
halos at such early epochs formed in C-I-HDM models.) 
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Fig. 7.— Redshift evolution of the physical correlation length vq for the simulated halos. In the top 
panels, the three curves from top down are for the = 0.3, 0.5, and 1.0 LCDM models. In the 
bottom panels, the two curves are for the = 0.2 (top) and 0.1 (bottom) C+HDM models. The 
redshift dependence of tq varies with cosmological parameters and is not a power law (see eq. [3]) 
as has often been assumed. The left and right panels compare two different halo mass ranges and 
show that more massive halos generally have a larger tq. 
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Fig. 8 .— The evolution parameter rj{z), defined as ro{z) = {l + z)~'^^^'^ ^o(O), for the M > Mq 

virialized halos in four models: CDM, flm = 0.5 LCDM, = 0.3 LCDM, and = 0.1 C+HDM 
models. The COBE normalization is used, including the CDM. The vertical labels on the right 
mark the corresponding e used in the literature, where slope of 7 = 1.8 has been assumed to relate 
r] and e: e = I .87 — 3 (see eq. [7]). The rise of 7 with decreasing z shows that the halo correlation 
length increases faster than a simple power law assumed in the literature. The clustering passes 
7 = 1 (i.e. e = —1.2) for fixed clustering in comoving coordinates and approaches 7 = 1.7 (e = 0) 
for stable clustering at z w 0 . 





